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Abstract 

The problem investigated in this paper is einselection, i. e. the selection of mu- 
tually exclusive quantum states with definite probabilities through decoherence. Its 
study is based on a theory of decoherence resulting from the projection method in the 
quantum theory of irreversible processes, which is general enough for giving reliable 
predictions. This approach leads to a definition (or redefinition) of the coupling with 
the environment involving only fluctuations. The range of application of perturbation 
calculus is then wide, resulting in a rather general master equation. 

Two distinct cases of decoherence are then found: (i) A "degenerate" case (already 
encountered with solvable models) where decoherence amounts essentially to approx- 
imate diagonalization; (ii) A general case where the einselected states are essentially 
classical. They are mixed states. Their density operators are proportional to microlo- 
cal projection operators (or "quasi projectors") which were previously introduced in 
the quantum expression of classical properties. 

It is found at various places that the main limitation in our understanding of de- 
coherence is the lack of a systematic method for constructing collective observables. 
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1 Introduction 



The discovery of decoherence has already much improved our understanding of quantum 
mechanics. The effect has now been observed experimentally pp. Many of its consequences 
have been obtained theoretically, but its foundation, the range of its validity and its full 
meaning are still rather obscure. This is due most probably to the fact that it deals with 
deep aspects of physics, not yet fully investigated. 

The intuitive idea of decoherence is rather clear The wave function of a macroscopic 
system depends on a very large number of variables and its local phase is very sensitive 
to boundary conditions, couplings and initial conditions. Any phase coherence between 
different components of the wave function is therefore exposed to destruction, after which 
macroscopic interferences disappear. It is unfortunately very difficult to build up a satisfac- 
tory theory on this intuition, because a knowledge of phases remains out of reach for the 
A^-body methods at our disposal. 

Some questions about decoherence 

The problems of decoherence are most often stated after making a few simple and prag- 
matic assumptions. One assumes particularly that a few collective (or relevant) observables 
can describe the main features of a (generally macroscopic) system, and they are known a 
priori. The system is then split formally into two subsystems: a "collective" one (which 
is associated with the relevant observables) and an environment, which can be external or 
internal. Each of these two abstract subsystems has its own Hilbert space. Tic or Tie and 
the Hilbert space of the whole system is the product Ti. = Tic ® "He- The Hamiltonian is 
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accordingly split into three parts, one for each subsystem and one for their coupling: 



Observers are supposed to have only a direct knowledge of the collective subsystem. Its 
properties are expressed by a reduced density operator p,,, which is obtained from the full 
density operator p through a partial trace on the environment 



The time evolution of pr exhibits the various aspects of decoherence. It has been investigated 
mostly by means of more or less exactly solvable models. Two models were particularly 
important because they were rather close to reality, at least in specific circumstances. In one 
of them the environment is replaced by a collection of harmonic oscillators [3-8]. Another 
model represents decoherence as an accumulation of scattering phase shifts when particles 
from an external atmosphere collide with a macroscopic object [H]. Much of what is known 
about decoherence was learned from these models and some of their variants jlUj . 

The conclusions have been accurately summarized by Zurek. He distinguishes three dif- 
ferent physical effects resulting from decoherence [lOj: There is first a destruction of macro- 
scopic interferences, then some privileged state vectors become selected as alternative phys- 
ical events and finally these states evolve classically. The privileged states are also called 
pointer states in analogy with the position of a pointer on a dial in a measuring apparatus 
Most models predict that these states exist and are orthogonal so that they define a 
Hilbert space basis in which the reduced density operator becomes approximately diagonal 
after a short while. The existence of this basis is essential since it defines a unique set of al- 
ternative events with well-defined probabilities. The name of "einselection" has been coined 
by Zurek for the mechanism selecting this basis. 



H = He® h + Ic® He + Hi . 



(1.1) 



pr = Tr^p . 
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These results have so far-reaching consequences for the interpretation of quantum me- 
chanics and other apphcations such as quantum computing that one must assess their exact 
degree of generahty. How much of them is specific to the models that were used and how 
much is universal ? This question raises several problems: 

• 1. A basic preliminary problem is concerned with the meaning of collective observables. 
When an actual physical system is given in practice, it is a rather straightforward matter 
to guess what coordinates describe most conveniently its macroscopic dynamics (the choice 
of these "generalized coordinates" goes back to Lagrange). But the question of defining 
correctly the collective observables for an arbitrary quantum state of the whole system, i. e. 
to select what is collective and what can be considered as an environment, is much deeper. 
It will be seen again and again in this paper that it represents the real limit of our under- 
standing. 

• 2. One may be puzzled by the fact that explicit models yield einselection somewhat too 
easily. This is because most of them rely on a coupling satisfying the commutation property 



between the coupling Hamiltonian Hi and a collective coordinate observable X (which may 
be supposed unique for simplicity). It is then found that pr becomes approximately diagonal 
in the basis |a; > of eigenvectors of X. It is clear however than Eq. ()1.3|) is very restric- 
tive, at least from a mathematical standpoint, and one cannot assume it to hold universally. 
What happens then when this condition is not satisfied ? Is there still some sort of diagonal- 
ization ? If so, along which basis ? To answer this problem will be the main task of this paper. 
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• 3. Problem 2 is made somewhat tricky because there exists a very large class of systems with 
collective coordinates for which condition ()1.3p holds. They are mechanical non-relativistic 
systems (excluding macroscopic electromagnetic effects), described by the position coordi- 
nates of hydrodynamics [12-15]. These systems, which are exceptional from a mathematical 
standpoint, may very well turn to be universal in measurement theory since a measuring 
device involves practically always some mechanical parts entangled with the rest of the 
apparatus. As will be shown in Section 7, the property ()1.3|) results from the form of ki- 
netic energies and two-body potentials in non-relativistic physics. This remarkable feature 
"explains" why classical mechanics can be formulated in ordinary three-dimensional space 
although the wave functions are defined on a configuration space [IT^ . 

So frequent an occurrence of a very special case may be puzzling from an intuitive stand- 
point. It raises a conflict between what we consider as general, either when speaking of the 
physical world or of the mathematics of the theory. This possible source of confusion will be 
avoided here by referring explicitly to "mathematical generality" when a "general property" 
or a "general case" will be mentioned, except when otherwise explicitly stated. 

• 4. One might be lured by models into premature conclusions and a sufficiently wide-ranging 
theory of decoherence is necessary for assessing general properties. The construction of such 
a theory is the fourth problem to be considered here. 

• 5. Finally, one must consider the attractive approach to einselection by Zurek [TUl I17j. 
Einselected states are supposed to be the most predictable (or robust) carriers of informa- 
tion. Given a collective state \Ef (which may be pure or not) and the corresponding initial 
density operator p*(0) = |\1' >< \E'|, one considers the time-dependent reduced density op- 
erator p^(t). Its ability to preserve an information content is characterized by some relevant 
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functional of which may be minus the von Neumann entropy or more conveniently the 
measure of purity 



This quantity is then used to construct a "predictability sieve" distinguishing among the 
states: The largest the quantity ()1.4|1 is, the more predictable the state \1/ is supposed to 
be. Model examples suggest that einselected states do exist and are rather insensitive to 
a change in the coupling or a redefinition of the environment. A fifth problem consists in 
evaluating this conjecture in a wider framework. 

The present results 

These five problems will not of course be solved here completely, but some definite or 
suggestive answers will be obtained. The most precise results are concerned with einselection 
and diagonalization, their meaning and their relation. As a preliminary, one needs a suffi- 
ciently wide-ranging theory of decoherence, as stated in Problem 4. The theory to be used 
here does not claim to be new. It relies on the familiar idea that decoherence is a special 
kind of irreversible process. This means that one may expect that the most general theory 
of decoherence presently at our disposal would derive from the existing quantum theories of 
irreversible processes. Moreover, the most convenient such theory is the so-called projection 
method [18-21]. Its main features are recalled in Section 2 and it is applied to decoherence 
in Section 3. Although this method was previously introduced elsewhere |15| 1^ 1^. some 
improvements will be required before applying it for the present purpose. These develop- 
ments are mostly given in the Appendices and they may be useful for using the method in 
other problems. 

A very important remark concerning this approach is the possibility of making a definite 



c^(t) = Trp%{t) . 
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choice for the couphng Hamiltonian Hi from which one can easily derive a master equation for 
Pr- The point is that, given a priori a couphng Hamiltonian, one can construct another (time- 
dependent) one consisting only of fluctuations. Standard perturbation methods can then be 
applied confidently in most cases. Although this procedure is familiar near equilibrium (it is 
used for instance in the derivation of the fluctuation-dissipation theorem it works also 
far from equilibrium as will be shown in Section 4, providing a strong handle on decoherence. 

The master equation one obtains in this way is probably the most general one that is 
accessible with present techniques and therefore the most appropriate one for investigating 
einselection, as done in Section 5. One thus finds that, contrary to current expectations, 
two different cases of einselection exist. The first one was encountered with solvable models 
and it is well known, although it is far from being general in a mathematical sense (in the 
case of an arbitrary coupling). It must be considered on the contrary as a degenerate case 
for the following reason: If n is the number of X-observables, decoherence is controlled 
in the general case by a differential Laplacian-like operator in the 2n- dimensional phase 
space. In the simplest case n = 1, the decoherence "Laplacian" is associated as usual with 
a quadratic form (like the two-dimensional Laplacian d"^ /dx^ + d"^ /dp^ is associated after a 
Fourier transform {x,p) —>■ {C,,i]) to the quadratic form (■C^ + z/^). The decoherence Laplacian 
is degenerate when it acts on only one variable (for instance the coordinate x and not the 
associated momentum p) so that the corresponding quadratic form is degenerate (having 
a zero eigenvalue). In the degenerate case, to which the usual models belong, decoherence 
essentially amounts to a diagonalization (in the basis einselected by degeneracy). 

The non-degenerate case is investigated in Section 6. The results are not those expected 
from Zurek's predictability sieve, at least as far as I understand it. There are generally two 
distinct times for decoherence. Typically, in conditions when these times are very different, 
decoherence selects a basis of privileged states in which pr begins to become diagonal, but 
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then these "einselected" states do not preserve their probabihties and they begin to share 
them with neighboring states. When the two decoherence times are similar, not much 
remains apparently of the idea of einselection. The final outcome of decoherence is rather a 
tendency towards uniformity where py. becomes as close to the identify matrix as its finite 
trace can allow (more precisely, it corresponds to a uniform Wigner function over a rather 
large region of phase space). Nevertheless, macroscopic interferences are still destroyed and 
classical behavior may follow. 

The question of the different time scales is discussed in section 7 and a strong connection 
is found with the special properties of hydrodynamical variables. The relation between 
decoherence and dissipation coefficients plays an important role in this discussion. Problem 
3 becomes then central because decoherence depends most often in practice on the fact 
that the space coordinates of a non-relativistic piece of matter satisfy the condition ()1.4|) 
implying degeneracy. One is thus led again to Gell-Mann and Hartle's ideas concerning the 
link between coarse graining and the existence of a diagonalization basis IT!?j . 

Finally, the occurrence of classical dynamics after decoherence is considered in Section 8 
in both the degenerate and the non-degenerate cases. In the non-degenerate case, although 
one can still speak of einselected states, they are far from being pure states. They are mixed 
states whose density operator at a definite time is proportional to a "microlocal projection 
operator", which is known in mathematics as the best expression of a classical property 
involving position and momentum together [25-27]. Finally, the conclusion goes back to the 
strong connection between the origin of classicality and Problem 1, i. e. the construction of 
collective coordinates. Some proposals for further research are indicated. 

Appendix A gives a derivation of decoherence theory from the projection method in 
irreversible processes. Appendix B is concerned with the relation between decoherence and 
dissipation coefficients. Appendix C shows how to include the insightful scattering model of 
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decoherence by Joos and Zeh in the framework of the present theory. 

2 A quantum theory of irreversible processes 

One needs a sufficiently wide-ranging theory of decoherence for asserting its general 
features. Since the loss of phase coherence through decoherence produces disorder, typical 
of an irreversible process, the most general relevant theory is the projection method, which 
is now briefly reviewed [18-21]. 

One considers a system with many degrees of freedom, whose density operator evolves 
according to the von Neumann-Schrodinger equation, 

P^~[H,P\. (2.1) 

Some observables are supposed to be particularly relevant for a given problem and one 
wants only to know the time evolution of their mean values. These relevant observables are 
denoted by AK Their set may be finite or not, countable or not. One usually includes the 
identity operator / in the set and also the conserved quantities such as the total energy H, 
although this recipe is not imperative. The "exact" average values of these observables are 

a^{t)=Tr(A^p{t)) . (2.2) 

The first step of the method consists in introducing a time-dependent test density oper- 
ator po(^) satisfying the following two conditions: (i) It gives the exact average values of the 
relevant observables: 

Tr {A^po{t)) = Tr [a^ p{t)) = a^{t) . (2.3) 

(ii) Its information content is minimal (which means that it does not provide anything else 
than the quantities {a^{t)}). It can therefore be written as 

Po{t)=exp[-X,{t)A^) (2.4) 
9 



where the numbers {Xj} are Lagrange parameters and the Einstein summation convention 
over repeated indices has been used. Since the identity / belongs to the set {A^, the trace 
of po is normahzed. 

One will use the name "density" in the present paper for an operator with a finite trace 
(also called a trace-class operator). It is neither supposed to have a unit trace nor to be 
necessarily positive. One defines a set of (time-dependent) densities 

Sj = dpo/da' . (2.5) 

They satisfy the important orthogonality properties 

Tr {s,A^) = 6{ , (2.6) 

amounting essentially to da^ /da^ = 61 in view of Eq. ()2.3|) . 

The theory makes use of "superoperators" , which act linearly on a density to yield another 
density. For instance, Eq. ()2.1|) can be written conventionally as 

p = Lp, (2.7) 

where L is the Liouville superoperator. Another important superoperator is defined by 

P = Sj(g) , (2.8) 

which means that when acting on a density p, it gives (with the summation convention) 

P/i = Sj ■ Tr . (2.9) 

It is a projection in so far as it satisfies the simple equation 

= P , (2.10) 

resulting from the orthogonality properties ()2.6|) . 
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One can define a "relevant" (time-dependent) density operator pi by 

Pi = Pp. (2.11) 

It generates the exact average values {a{t)} since 

Tr (aVi) = Tr [A^Pp) = Tr (^A^Si) Tr [A'p) = 6{a^ = a? . (2.12) 

(As a matter of fact, it coincides with the test density operator po)- 

Denoting by X the identity superoperator, one also introduces Q = I — P , which satisfies 
the projection property = Q in view of Eq. ()2.10|) . as well as the orthogonality properties 
QP = PQ = 0. One defines then another density p2 = Qp (so that p = pi + p2 )• Applying 
the two projections P and Q to the evolution equation ()2.7|) and taking into account the 
time dependence of these projections, one obtains two coupled evolution equations 

pi = PLPpi + PPpi + PLQp2 + PQp2 , (2.13) 

P2 = QLQp2 + PQp2 + QLPpi - PPpi . (2.14) 

A last step would be to eliminate p2 to obtain a master equation for pi but it will be left for 
the special case of decoherence. 

3 The case of decoherence 

The previous theory can be now applied to decoherence and some preliminary consid- 
erations will make the task clearer. A first problem is to choose the relevant observables. 
If one thinks of macroscopic interferences, it is clear that they can involve many different 
macroscopic observables so that every collective observable is relevant. When the collective 
observables describe a measuring apparatus, the measured microscopic observable is also 
relevant although not collective. 
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The environment can be defined by its observables, which commute with the collective 
ones. The resulting splitting of the system into a collective one and an environment is 
generally time dependent since for instance every new bubble in a bubble chamber or every 
new spark in a spark chamber generates new collective observables. Decoherence is such 
a rapid process however that the collective environment splitting can often be considered 
as fixed during the very short time of an individual decoherence process, justifying the 
expression (jl.lj) for the Hamiltonian. The coupling Hi in this equation is responsible for the 
interactions between the collective system and the environment, including dissipation and 
decoherence. 

The set of relevant observables is completed by the identity operator / (insuring nor- 
malization) and the environment Hamiltonian H^, or more properly 1^® He- The total 
Hamiltonian H might have been used as relevant in place of but this choice would have 
been inconvenient, as will be seen later). 

Introducing an arbitrary orthonormal basis |A; > in the collective Hilbert space, the set 
{\k >< k'\} provides a linear basis for the collective operators. A look at the calculations in 
Section 2 shows that they nowhere use the fact that the {A^} are hermitian so that one may 
use the set of operators {\k >< k'\} as relevant "observables". Alternatively, one might use 
the set of hermitian operators 

{(1/2)(|A; > ±\k' >)(< k\± < k'\); (1/2)(|A; > ±t\k' >)(< k\Ti < k'\)} 

as relevant with identical results. Anyway, the set of relevant observables for a theory of 
decoherence will be chosen as 

1^'='^' = \k><k'\® Je, = I^A^ = h® He) . (3.1) 

Since none of them connects the collective and the environment Hilbert spaces, the test 
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density operator ()2.4p is a tensor product 

P0 = Pc® Pe ■ (3.2) 

Applying Eq. ()2.3|) to the operators A''^'^ one finds that 

< k'\pc\k >= Tr [a'^'^'po) = Tr (a'^'^'p) =< k'\trp\k >=< k'\pr\k > (3.3) 

so that the collective test density pe is identical with the conventional reduced density p.2|) . 
The second equahty results from Eq. ()2.1|) . A convention for traces that will be used 
everywhere has also been introduced, the notation Tr standing for a trace on the full Hilbert 
space and tr for a trace on the environment. 

According to Eq. ()2.4|1 . the environment test density pe is given by 

Pe = exp {-PHe - a) , (3.4) 

where a is a Lagrange parameter insuring normalization and (3 insures that the energy 
Hf. of the environment has its true average value E. This density is the same as if the 
environment were in thermal equilibrium but it should be stressed that it is only an auxiliary 
mathematical quantity providing a correct (time-dependent) value for E with no assumption 
about equilibrium. 

In Appendix A, the auxiliary densities Sj are obtained from Eq. ()2.5j) . Denoting respec- 
tively by Si and S2 the densities associated with / and H^, one gets 

Skk' = \k' >< k\(g) Pe , (3.5) 

Se = Pc® Pe {He - E) A^^ ^ (3,6) 
Si = -Epe ® Pe {He ~ E) A'^ , (3.7) 

where A is the uncertainty in energy 

A^ = tr (Hlpe) - E^ . (3.8) 
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When acting on an arbitrary density /i, the projection P is given according to Eqs. ()2.9|) 
and (|2.nil2.7|l by 

Pfi = trfi ®p, + (pe ® {pe{He " E)A-^}) ■ {TrH.fi - ETrfi) , (3.9) 
from which the relations = P and pi = po follow. 

Specifying the coupling 

One may now introduce an important remark that will later justify the use of pertur- 
bation theory. To begin with, one may notice some arbitrariness in the splitting of the full 
Hamiltonian H into three different terms as in Eq. (jl.lll . A simple recipe for fixing them is 
to impose that 

trH^Pe = . (3.10) 

The meaning of this condition can be seen on the example of a cylinder containing a gas. 
A collective coordinate is specified by the position x of a piston whereas the environment 
consists of the gas and the matter of the piston itself. A straightforward definition of Hi 
could be the sum of the potential energies between the atoms in the piston and the gas. This 
interaction is far from being weak, since the confinement of a gas is not a weak effect, but a 
large part of it consists of a collective energy since trHiPe is a collective operator. One can 
then change the definition of the different parts in H by including this operator in He and 
removing it from Hi or, more precisely, by introducing 

H'^ = H, + tr{Hip,) , 
H[ = Hi- tr {Hipe) ® /e . 

The quantity trHip^ represents in this example the effect of the gas pressure on the piston. 
The new expression of H^. is time-dependent (like pressure) but the new expression of Hi 
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satisfies the condition ()3.10|) . It consists only of tlie pressure fluctuations resulting from the 
collisions of the gas molecules with the piston. 

The idea of introducing a purely fluctuating coupling and to use perturbation theory for 
computing its effects is familiar in quantum fluctuation theory [21]. The fact that one can 
still use it far from equilibrium when dealing with decoherence is due to the possibility of 
representing everything collective by the test density. From there on, the condition (|3.1Up 
will be assumed. 

One may also understand at this point why the choice of He as a relevant observable is 
more convenient than the total Hamiltonian if, which is usually recommended This 
is because the expression ()3.2|) for the test density implies the simple rule fl3.10|) for the 
coupling, with the benefits just mentioned. Everything would have been more obscure and 
would have implied much heavier calculations if H had been chosen as a relevant observable. 

The evolution equations 

It is easy to write down explicitly the evolution equations ()2. 13112. ll)) for the case of 
decoherence. It is convenient to split Eq. ()2.13j) for pi into an equation for and another 
for pe (or for the time evolution of the internal energy). This is done by taking respectively 
the traces of Eq. ()2.13|1 on the environment and the collective Hilbert spaces. As shown in 
Appendix A, the results are: 

Pr = -'-{[He,pr]+tr[H,,p2]) , (3.11) 

E+'-Tr{H,[Hi,pi + p2])=Q . (3.12) 
As for the second evolution equation ()2.14|) . it becomes 

P2 = -(V^) [H, Pi + P2] + ii/h)tr i[Hi,p2]) ®Pe-pr®Pe. (3.13) 
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4 A master equation 

The most delicate step in the projection method consists always in "solving" the second 
evolution equation ()2.14|) for p2 in terms of pi before inserting the result into Eq. ()2.13|) 
pij . This is much easier when perturbation theory can be used. Perturbation calculus 
has been used already in the present framework when Hi is known a priori to be small, 
as often happens in quantum optics [221 12^1 • It should also presumably be valid in many 
instances when condition ()3.10|) is applied and Hi is a pure fluctuation (although one must 
acknowledge that a purely fluctuating coupling does not insure with certainty the validity of 
perturbation calculus). Anyway, according to Appendix A, the evolution equations (3.11-13) 
become at leading order in Hi : 

Pr = ~ {[H,,p,]+tr[Hi,p2]) , (4.1) 

P2 = -{i/n) [Ho, P2] - m) [Hi, P2] . (4.2) 
In the second equation, Hq denotes the uncoupled Hamiltonian 

HQ = H,®h + h® H, (4.3) 

and Eq. ()4.ip is exact whereas Eq. ()4.2|) is valid at flrst order in perturbation theory. 
The second equation is easily solved after introducing the evolution operator 

U{t) = exp {-iHot/h) . (4.4) 

Strictly speaking, Hq is generally time dependent and the integrand in Eq. ()4.4p should be 
replaced by an integral of Ho{t) on time. It is difficult however to conceive of a case where 
this external time dependence is not much slower than decoherence and the expression (j4.4p 
is therefore most often valid as it stands. If not, the necessary changes are so trivial that 
they need not be mentioned here. One thus gets 
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P2{t) = -{i/h) f dt' Uit-t')[H,,p,{t')]U-\t-t') 



(4.5) 



J —oo 



No effect of the initial value of p2 (at time — oo) has been included in Eq. ()4.5j] . This is 
justified when the environment is initially in thermal equilibrium (since then p2{—oo) = 0). 
More generally however, it may be expected that an initial lack of equilibrium does not 
influence the decoherence effect, so that Eq. ()4.5|1 is valid for our present purpose. This 
point was checked in a special case by Paz, Zurek and coworkers |28|, 1^ . 

Inserting Eq. ()4.5j) into Eq. ()4.H1 . one obtains the following "master equation" for 
decoherence 



The first term in the right-hand side represents the quantum evolution of the reduced density 
operator under the action of the collective Hamiltonian He- The second term is responsible 
for decoherence. 

This equation is not new but it was derived previously either under the assumption of 
a small coupling |22l 12^] , or as a guess P^J . It will be used here as a sufficiently general 
framework for a study of einselection. 

The wide range of this master equation is confirmed by its agreement with previous 
models. This is easily shown when the environment is represented by a collection of harmonic 
oscillators [3-8]. The key experiment by Brune et al. showing the existence of decoherence 
also belongs to the domain of Eq. ()4.6p since Hi is small in that case ED] • In the case 
of the collision model by Joos and Zeh the calculations are less trivial and they are given in 
Appendix C as a non-trivial example of the master equation universality. 




dt' tr Hi, U{t - t') [Hi, pr{t') ® peit')] U-\t - t') 



(4.6) 
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5 Decoherence versus diagonalization 

Models have been extremely useful for understanding decoherence. When the collective 
subsystem is described by a few position observables X, decoherence was found to diagonalize 
the reduced density pr, in the basis \x > consisting of eigenvectors of X. The question to be 
now considered is therefore: Does decoherence always implies some sort of diagonalizationn 
? Is there always a selection of privileged "pointer states", or einselection as defined by 
Zurek ^ ? 

One may first select the equation on which this question will be investigated. The idea 
of diagonalization must be used with some care because the reduced density operator never 
becomes completely diagonal in view of the first term in Eq. ()4.6p representing collective 
dynamics. For a finite value of the difference x — x' the matrix elements 

Pr{x,x';t) =< x\pr(t)\x' > , (5.1) 

vanish exponentially with time, whereas microscopic values of x — x' are dominated by 
collective dynamics and they remain finite. This is why there is decoherence on large scale 
while atomic physics remains perfectly valid at small scale. The question of diagonalization 
is therefore much clearer if one leaves aside the first term in Eq. ()4.6p and consider "pure 
decoherence" as the behavior of a density operator obeying the truncated equation 

p^ = -(^l/n'^) t dt' tr\Hi,U{t-t')[Hi,pr{t')®pe{t')]U-\t-t')]=D . (5.2) 

The main task will then consist in an analysis of the right-hand side of Eq. ()5.2|) ). which 
has been denoted by Z^. It will also be useful to introduce the notation 

Hi = U{t - t')HiU-\t - t') , (5.3) 

so that one has 

D = f dt' tr \Hi, [iff, Uc{t - t')pr{t')U-\t - t') ® pe{t')\\ (5.4) 

J —oo '- '- 
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(where ^ PeUe has been replaced by pe in view of Eq. ()3.4|) . 



Weyl symbols 

One will consider the case when there exists a set of n collective "position" observables, 
altogether denoted by X. The quantity D is itself a collective operator and it will be 
convenient to describe it by means of a Weyl symbol , in analogy with the description 

of pr by a Wigner function [S2]- The standard Weyl calculus can be slightly generalized to 
include "operator-valued symbols" acting on the environment as follows: 

Let A denote an arbitrary operator in the full Hilbert space (such as Hi for instance). 
Introducing the basis {|x >} in the collective Hilbert space and an orthonormal basis {|^ >} 
in the environment Hilbert space, the matrix elements of A can be expressed through a 
partial Fourier transform 

I X -\- x' \ 

{27ch)~"'dpAnn' I — - — , p j exp {ip{x' — x)h} . (5.5) 

Every quantity (x, p) is a function of (x, p) and the ordinary Weyl-symbol of the matrix 
element Ann' =< n\A\n' >, which is a collective operator. It will be convenient to consider 
it as the {n, n') matrix element of an operator- valued symbol A{x,p), which is a function of 
{x,p) and an operator in the Hilbert space of the environment. 

The symbol of the product AB of two operators A and B can then be expressed as a 
series in powers of h involving their symbols jSH 123 '■ 

AB = A-B- {ih/2) (ApB, - A,Bp) - {h^ /2A) {a^^B,^ + A,2Bp2 - 2Ap,Bp,) + Oih') . 

(5.6) 

The notation has been simplified by omitting the arguments {x,p) of the symbols and lower 
indices stand for derivatives (for instance, A^p stands for d'^A{x,p)/dxdp). Eq. ()5.6|) is well 

19 



known in Weyl's calculus when the symbols are ordinary functions. It is easily extended 
to operator-valued symbols by considering matrix elements and a unique new rule must 
be added to the case of functions: the order of the operators in the product AB must be 
respected in the products of symbols and their derivatives. The symbol of the reduced 
density operator p,. is the Wigner function, which will be denoted by W{x,p). It is not 
operator-valued and commutes with other symbols. 

The only further formula one will need from the Weyl calculus is the expression of a 
complete trace: 



Calculation of the decoherence term D 

It will be convenient to consider from there on the case of a unique coordinate X {n = 1) 
although the generalization to arbitrary values of n is trivial. Applying Eq. ()5.fi|) to the 
double commutator in Eq. ()5.4|1 . one obtains the symbol D of the decoherence term D at 




(5.7) 



order , as shown in Appendix B: 



D 



dt' (d/dx (C^^VTj + C^^Wj) + 3/ dp [C^^'W^ + Cf^Wj)) . (5.8) 



The function W 



{x,p) is the symbol of the collective operator 



Ucit-t')prit')U;\t-t') . 



(5.9) 



The various decoherence coefficients are given by 




(5.10) 



(5.11) 



(5.12) 



(5.13) 
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It is possible in principle to derive the main consequences of the master equation for 
decoherence from these equations by using the powerful methods of microlocal analysis 
PH] . It will be much simpler however to rely on a few usual approximations. The first one 
assumes that the coefficients (5.10-13) vary slowly with {x,p) or, more precisely, one neglects 
the collective evolution Uc{t — t') in the factors U and occurring in the expression ()5.3|) 
of Hj^. The physical meaning of this approximation is discussed in Appendix B, where the 
following expression of C^^ is obtained: 

(JPP ^ J2 Hi^nn' Hi^n'n exp {iuJnn'{t - t'))pnn' COsh {pkuJnn'/^) , (5.14) 
nn' 

where the states \n > are the energy eigenstates of Hf. with eigenvalues £"„ and one has 
written 

Hlxnn' =< n\dHi{x,p)/dx\n' >,UJnn' = {En - E^') /h,Pnn' = exp {En + Cn') 1'^,- a] . 

(5.15) 

Eq. ()5.14p suggests that the relevant frequencies ujnn' in the sum are contained in an 
interval [— i^, Vt] characterizing the environment and generally large as compared with the 
rate of collective dynamics {Vl is typically a Debye frequency for an internal environment). 
Hu, Paz and Zhang have shown that the master equation is instantaneous (i. e. involves 
no retardation) in the case an oscillator environment, when the collective Hamiltonian also 
describes an oscillator ^2j- This is due to the linear character of the equations in that case 
|33j . The resulting master equation has been solved explicitly by Ford and O'Connell jS]. 
This situation is however exceptional and the neglect of retardation is almost always an 
approximation. The question of its justification is tricky and it would warrant a separate 
investigation. When retardation effects are unimportant anyway, the integration on t' in 
Eq. ()5.8|1 is performed as if the integrand were a delta-function in time. The time-delayed 
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function is replaced by the ordinary Wigner function W and Eq. ()5.8p becomes 

D = ^ {g^^W, + g'^'Wp) + ^ (gPW, + g^W,) . (5.16) 
The new coefficients are given by 

g'^ = f C'^{t-t')dt' (5.17) 



(with indices (i, j) = (x,p)). Exphcit expressions of these coefficients are given in Appendix 
B, showing that the coefficients g^^ and g^^ are positive symmetric : g^"^ = g^^, and the 
quadratic form 

c/^^a^ + 2g''Pa(5 + g^^p^ (5.18) 

is non-negative. One must then distinguish two significantly different cases according to 
whether the form ()5.18|) is degenerate or not, i. e. whether the determinant g^^gPP — {g'^^Y 
is zero or positive. 

The degenerate case 

The degenerate case was encountered in most models and only one coefficient, namely (7^^, 
was different from zero. It is then convenient to go back to the matrix elements pr{x,x';t) 
by inverting the Fourier transform ()5.5j) so that the pure decoherence master equation ()5.2|1 
becomes 

— Prix, x'; t) = ■ " x'Y pr{x, x\ t) . (5.19) 

Diagonalization in the basis {|x >} is then obvious when g^P is a constant since the solution 
of this equation is 



Pr{x, x', t) = Pr{x, x' , 0) CXp 



2 - ^') t 
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Similarly, when the only non-zero coefficient is g^^, one may use the momentum basis {|p >} 
to obtain: 

— p,.(p,p';t) = — ^ ■ {p-pfpr{p,p']t) , (5.20) 

implying again diagonalization. 

A simple condition for the coupling implying diagonalization in the position basis is 
given by Eq. (jl.Hj) |311 1^ . Using coarse graining, Gell-Mann and Hartle have shown that 
this condition is satisfied for mechanical systems when using hydrodynamical observables as 
relevant [T!^ . 

6 The non-degenerate case 

Quite different results are obtained in the general case when there is no degeneracy. One 
may note first that the differential operator in the right-hand side of Eq. is similar to 

a Laplacian, which is given by 

in the case of a metric ds^ = Qijdx^dx^ (with Qijg^'' = 6f). The factor g is the determinant 
of the matrix with elements Qij or the inverse of det ((?*•'). One could use this remark in 
principle for a general study of decoherence but it would need the full power of microlocal 
analysis. Rather than entering into such heavy mathematics, it will be convenient to consider 
only the case when the coefficients g^^ are constants. A further simplification is obtained 
by diagonalizing the quadratic form (jS.lSj) . This is done by a change of variables after 
introducing scale-invariant parameters: Let L be a unit of "length" (i. e. a scale with the 
dimensionality of X) and 11 a unit of momentum. The transformation 

UX' = UX cos e + LP sin 9, LP' = -UX sin 6 + LPcoO , (6.2) 
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can be viewed either as an "orthogonal" change of axes in the {x,p) plane or as a linear 
canonical transformation. It leaves Weyl's calculus invariant so that if one chooses 6 to 
diagonalize the metric, one obtains a simpler equation for pure decoherence, namely (after 
dropping the prime indices) 

dW d^W d'^W 

General decoherence is not a diagonalization process 

The general case of decoherence occurs when the quadratic form ()5.18|1 is non-degenerate. 
Does then the effect still amount to diagonalization ? By looking at the degenerate case, one 
sees that diagonalization was due to a specific property of the collective operator D: There 
was a specific orthonormal ("pointer") basis {|j >}, such that 

< j\D\j >= for each j , (6.4) 

Re < j\D\k >< , for every pair of indices j ^ k . (6-5) 

These relations held true for any density matrix pr entering in the definition of D. They 
must obviously be satisfied if diagonalization takes place, at least if the basis is independent 
of the preparation Pr(0) and depends only on the decoherence coefficients. They do not hold 
however in general as shown by the following 

No-go theorem 

Whatever the state ip, it is impossible for the equation 

< ilj\D{pr)\^p >= (6.6) 

to hold for every density matrix pr, 
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Proof According to Eq. ()5.16p , one can write 

D{p) = Ap , (6.7) 
where A is understood as a superoperator acting on a collective density p. One can write 

<ip\Ap\il: >=Tr^{\ilj ><ij\Ap) . 

If this equation is supposed to be valid for any choice of p, one must have (since the super- 
operator A is hermitian) 

A\iIj >< = . 

When written explicitly in the position basis, this equation becomes 
r 02 



from which one gets 



^(x + ^/2)r{x - m = a(0 exp (^^^{x^/h)j + 6(0 exp (^-^^{x^/h)j . 

This is however impossible (even if the coefficients are distributions) because it would imply 
that the wave function of the state ip increases exponentially for large values of its argument. 

Note : The present theorem forbids the existence of a universal diagonalization basis. The 
possibility of a p-dependent basis remains open, although it looks very doubtful. 



Decoherence in the non- degenerate case 



Since decoherence cannot be generally a diagonalization process, one must investigate 
it anew. Its consequences are most easily obtained when x — x' is large. It will be more 
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convenient to use the notation {x',x") for the arguments of the reduced density matrix 
Pr{x',x") in the position representation and to introduce auxihary variables x = {x' + x")/2, 
^ = x' — x" . This means that we are interested on the case where ^ is large (macroscopic). 
After performing a Fourier transform to go back from the variable p to S,-, the pure decoherence 
equation ()5.2|) becomes 

p = 9^'^ - {re/h') p . (6.8) 

The time evolution of the function p(x, ^) = pr{x + ^/2,x — C,/2) is therefore given by 

p{x',^,0). (6.9) 

The first factor on the right-hand side shows that pr{x + C,/2,x — ^/2) tends to become 
diagonal in the position basis, as in the degenerate case. The heat kernel in the integral has 
however a very different effect since it smoothes off the reduced density along the diagonal, 
so that probabilities that were initially distinct become mixed together. If the process is 
stopped at some time t, its effect is analogous to an imperfect measurement of the position. 

The smoothing effect is most clearly seen by considering as initial state a superposition 
of two distinct wave functions : 

p^(t = 0) = \ijj >< '01 , with l^jj >= lipi > +\^/j2 > (6.10) 

the two wave functions ipi{x) and ip2{x) being clearly separated with clearly different average 
values for X or P or both. One is interested in the interference part of pr originating from 
>< il!2\ and \iIj2 >< "^il in the initial state operator, but one must also now consider the 
probabilistic part originating from {ipi >< ipil and \ip2 >< 4'2\- The first factor in the right- 
hand side of Eq. ()6.9|) suppresses very rapidly the interference terms in p^ when ipi{x) and 
ip2{x) have well-separated mean values of the position X. It also suppresses them, although 
less rapidly when the values of < X > coincide while those of < P > are significantly 
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pix, t) = exp (-g^^et/h' 



X 



dx' exp -{x - x! f l^g^^'t 



different From tlie standpoint of macroscopic interferences, there is therefore nothing 
new. 

The smoothing integral in Eq. ()fi.l()|) introduces a new effect. It mixes together the 
probabihties for different values of < X >. If the state ()6.10p represents for instance the 
state of a pointer after a measurement, two results that would be distinct for an apparatus 
with degenerate decoherence can become indistinguishable if decoherence is non-degenerate. 
This conclusion does not depend on the specific form of ipii^) ^iid il'2{x)- It also holds for 
coherent states, which are not therefore properly einselected. 



A symmetric form of decoherence 



A convenient expression of non-degenerate decoherence can be obtained for any number 
n of collective variables. One denotes altogether by {^■'} the set of the 2n position and 
momentum variables and the pure decoherence equation becomes 

d ( ■■dW\ 

If the decoherence coefficients g'^^ are constants, one may introduce the inverse "covariant" 
coefficients gij satisfying the relations g^-'gjk = ^k- They exist only in the non-degenerate 
case. The solution of Eq. is then given by 

Wm = {ntr'J ^d^^r]exp[-g,k{e - V') - V') /^W{r],0) , (6.12) 

where g is the determinant of the matrix with coefficients gij, inverse of the matrix of the 
decoherence coefficients g^^ . It may be useful to notice that Eq. ()6.12j) remains approximately 
valid when the decoherence coefficients are not constants but slowly varying 

The effect of decoherence is therefore to smear out the Wigner function in phase space. In 
this approach, the removal of interference terms is due to the fact that an interference term, 
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localized in phase space with a vanishing integral, is rapidly reduced to zero under smoothing. 
One may also notice that the results are unchanged under a linear canonical transformation 
in phase space, at least when the coefficients are constants. They are only slightly modified 
when the coefficients are slowly varying under a smooth canonical transformation 

7 More about degeneracy 

One may now consider the order of magnitude of the decoherence coefficients. It should 
be stressed first that the condition ()3.1|) implying diagonalization was most often imposed a 
priori in the construction of models. If one again considers a unique position observable X, 
only one decoherence coefficient, g^^ , is different from zero when this condition is satisfied (as 
shown in Appendix B). Models have revealed a strong connection between the decoherence 
coefficient g^^ and the friction coefficient 7^^, which appears in the classical limit of the 
equation of motion when the collective Hamiltonian is = /2m + V{X). Classical 
motion is then governed by 



As shown in Appendix B, the spectral densities of the coefficients g^^ and 7^^ are very 
similar and, at high enough temperature, the two coefficients have a simple proportionality 
relation : 



What should be considered in that high temperature has been discussed by Hu, 



of decoherence implies a strong effect of decoherence as soon as 7^^ is not very small, i. e. 
when there is a possibility of dissipation. 



dp 



dV/dx - Y^p . 





Paz and Zhang [7j. The fact that g^^ enters with a denominator in the expression ()6.9|) 
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Rough orders of magnitude 



Let us consider a model with an environment of oscillators, the collective system being 
also an oscillator with frequency uj [3-8]. One assumes usually a coupling proportional to 
the collective position X: 



so that the condition p.3|) is satisfied. More generally, one may consider a coupling with the 
creation and annihilation operators of the collective oscillator in place of X, i. e. 

Hi = {X- iP/muj) ■ J2 + {X + iP/mu) ■ ^ \*a\ . 

i i 

According to Appendix B, one has then typically: 

According to Eq. the decoherence time characterizing the vanishing of non-diagonal 

interference terms separated by a distance Ax is given as usual by 

whereas the characteristic time after which there the probabilities are mixed up for two 
different positions on the diagonal separated by the same distance Ax is 

For reference, it may be recalled that the time necessary for the spreading of a wave packet 
on the same distance is 

TTlAx'^ 

twp = — ^ — • (7.6) 
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For not too small values of Ax, and a generic coupling (i. e. no degeneracy) one has for the 
various rates: decoherence ^ probability mix up ^ wave packet spreading. 



On the existence of pointer states 



The main conclusion of the previous sections was that the existence of an exact diago- 
nalization basis is not essential for most physical consequences of decoherence. On the other 
hand, it will be now shown that there is a very large class of physical systems for which such 
a basis exists. 

Coming back to the case when X denotes a class of n collective coordinates, there is 
exact diagonalization in the basis \x > if the derivatives of the symbol Hi with respect to 
the canonically conjugate variables p vanish. One has then according to Eqs. ()5.6|) and 
flB.9IIB.lljl : 

[Hi,X] = 0. (7.7) 

Under the same assumptions, according to Eq. ()B.17IIB.T9|) . the friction coefficients 7^^^ 
and 7^^ also vanish. 

When they do not vanish, the classical equations of motion become 

dxi/ dt = Pi/ rrii - 'jffpj - 7f/xj , (7.9) 
dpi/dt = F,--f- tf% - -ff^'xj , (7. 10) 

where Fi denotes a force. 

These equations look rather unusual and it is important to understand why they are 
exceptional (or unrealistic). In the case of a mechanical system (with no macroscopic elec- 
tromagnetic effects) we are familiar with a unique type of friction coefficient (7^^ ) and of 
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decoherence coefficient [g^^). The reason why was clearly shown by Gell-Mann and Hartle 
|13j . who used as position observables the hydrodynamical variables resulting from a coarse 
graining. The corresponding variables can be identified with the center-of-mass positions Xi 
of small pieces of matter, small enough from a macroscopic standpoint although containing a 
large number of atoms. The key feature is then the non-relativistic form of the Hamiltonian 
for the particles of matter, 

where the summations are performed over the particles (indicated by Greek indices) . If there 
is no macroscopic magnetic field (so that one can neglect the magnetic potential A), one of 
the Heisenberg equations of motion, yields the following simple relation between the classical 
velocity and momentum (denoted by Latin indices) 

Xi = Pi/rrii . 

Comparing this with Eq. ()7.9p . one sees that 7^^' = 7^^ = and, from Eq. ()B.16|) . one may 
expect that g^^ = g^^ = 0. Eqs. ()7. 9117. 10]) strongly suggest that this property follows from 
the Galilean invariance of non-relativistic mechanics under a change of reference system. 

Strangely enough, no realistic example of the non-degenerate case has yet been proposed, 
except tacitly in unprecise measurements. Examples might be expected however in electro- 
magnetic systems (where the magnetic and electric fields replace position and conjugate 
momentum) but the prospect of producing quantum superpositions of fields and see their 
decoherence seems rather remote. One must probably attribute the rarity of examples to 
the fact that decoherence has been mostly studied in the framework of measurement theory. 
There is almost always (or always) in that case some mechanical part of some apparatus 
that is entangled with the measurement result and the rest of the system, and it enforces its 
own einselection on them. 
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8 Classical behavior 

The derivation of classical behavior 



Decoherence is in most cases immediately followed by a classical behavior of the collective 
subsystem jH [TUl [121 1231 • Although this property will not be analyzed in detail in the present 
work, a few points involving again the problem of einselection are worth mentioning. 

Decoherence is described in the non-degenerate case by Eq. ()6.12|) involving a smearing 
effect on the Wigner function W{x,p). An important consequence is to make this function 
non-negative so that its interpretation as a density probability in phase space becomes signif- 
icant |MJ. As far as orders of magnitude are concerned, one may consider that a derivative 
operator d/dx acting on W is of the order of (g^^t)"^^"^ for t large enough (i. e. when 
decoherence is effective) whereas d/dp is of the order of (g^H)^^^'^. One may then consider 
more carefully the first term in the master equation ()4.6|) giving the following contribution 
to the master equation 

Pr = -^[H,,Pr] . (8.1) 

One can write down this equation in terms of the Wigner function and the Hamilton 
function h{x,p), which are respectively the symbols of pr and of He. It reads to second order 
in h : 

dW dhdW dhdW h"^ (d^hd^W d^h d^W d^h ff^W d^hd^W] 
dt dp dx dx dp 24 | dp^ dx^ dxdp^ dpdx"^ dpdx"^ dxdp"^ dx^ dp^ J 

(8.2) 

One recognizes in the first term in the right-hand side a Poisson bracket of the Hamilton 
function and the probability density, in agreement with classical physics. This term generates 
a classical evolution of the Wigner function as if its arguments (x, p) were moving according 
to the classical Hamilton equations. 
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The order of magnitude of the Ty-derivatives resulting from Eq. ()6.12|) imply that higher 
order corrections in h are negligible so that after some decoherence the evolution becomes 
classical. It is somewhat paradoxical that the analysis is more involved in the degenerate 
case PO]. The difficulty arises from a linear superposition of two initial wave functions with 
the same a;-location but different average values of p. The destruction of interferences must 
then wait till the motion due to the difference in the values of < P > separates the wave 
functions in x-space. A conspiracy of decoherence with the collective dynamics is therefore 
necessary for producing finally a classical behavior. 

Classicality and the choice of a collective subsystem 

The previous conclusion of a classical behavior assumed tacitly that the derivatives of the 
Hamilton functions are not large, but one might then get involved in a circular argument. The 
collective observables are chosen ordinarily on empirical grounds, from a direct knowledge of 
the system. One says: "I look at the system and I clearly see how it can be described by some 
coordinates, which I replace by quantum observables" . Then one concludes after much work: 
"See! The description of the system with these variables becomes finally classical" . This is 
certainly not a proof of classicality resting on the basic principles of quantum mechanics, but 
only a check of consistency: Classical behavior can be proved when the convenient variables 
for describing it have been selected by means of one's classical intuition. 

The question "How does one select a collective subsystem?" is therefore prior to the 
question of einselection. One can then look at Eq. ()8.2|1 from a different standpoint: It 
should give a criterion for choosing the collective observables and not provide a proof that 
they describe a classical motion. This criterion implies that the derivatives d/dx and d/dp of 
the classical Hamilton function in the second and higher terms of Eq. ()8.2|1 are not controlled 
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by factors involving some power of h~ . 

This gross criterion can be presumably much refined in view of Fefferman's formulation 
of quantum mechanics through pseudo-differential calculus (microlocal analysis) |SZj. He 
investigated the eigenstates of the complete hamiltonian H of an arbitrary quantum system 
by analyzing its symbol H in the phase space of the constituent particles and, by cutting 
this space into "distorted boxes", he was able to diagonalize H approximately. This is a 
deep result of abstract mathematics but there has been no direct application of it in physics. 
Nevertheless, it means that there exists one (or several) privileged ways of cutting phase 
space, into well-defined boxes, according to the possible states of the system. One could 
then envision that a pair of variables {x,p) is collective if it defines locally a 2-dimensional 
plane along which H varies very slowly. Such a property is strongly suggested by Feffermann's 
construction and it agrees with the small derivatives we just found characterizing classical 
behavior after decoherence. One may also presume the existence of a whole hierarchy of 
collective 2-directions, which would be ordered according to the magnitude of the derivatives. 

I will not try to elaborate further on this idea, which was proposed some years ago al- 
though not much progress has been made since j^H]- It represents however an alternative 
to Zurek's predictability sieves (with which it may be related). In any case, it stresses 
again that the most important problem in a real understanding of decoherence is an explicit 
construction of the collective observables (with a corresponding explicit definition of the en- 
vironment). This is closely related with a search for a real theory of the Heisenberg frontier, 
as also noticed by Zurek [TUj . 

Einselection 

Zurek's concept of predictability sieves was applied successfully to the case of an under- 
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damped collective oscillator interacting with an environment of oscillators [17,40-43]. It was 
found that Gaussian pure states with various average values {x,p) of position and momen- 
tum are selected in that case as the best carriers of information, suggesting more generally 
that some sort of coherent states would be einselected by decoherence just before classical 
motion. There is something puzzling however in the fact that the width of these Gaussian 
states is controlled by the parameters (m, u) of the collective oscillator (it has the same 
width as the ground state wave function of the oscillator). When looking at Eq. ()6.12j) . 
one finds on the contrary that decoherence in the non-degenerate case is insensitive to the 
characteristics of the collective Hamiltonian and it is completely determined by the coupling 
Hamiltonian through the decoherence coefficients. 

It may be recalled in this connection that a convenient family of einselected states was 
proposed earlier, although this name was not used [SHlEni- These states are closely related 
to Hormander's notion of microlocal projection operators. The symbol P{x,p) of such an 
operator is zero outside a regular cell C in the (x, p) phase space (i. e. a cell whose volume 
and boundary shape have large characteristic dimensions in terms of the Planck constant). 
P{x,p) is equal to 1 in C, except near the boundary where it goes smoothly to zero. The 
corresponding operator P is practically a projection [25] . 



Consider then the integrand of Eq. ()6.12j) for definite values of t , ^ = {x,p) and 
rj = {x',p') 



It can be considered as the symbol of a density operator Z originating at a time t large enough 
through non-degenerate decoherence from an initial state localized in the neighborhood of 
1]. Using Eq. ()5.6|) . one finds that PZ = Z if the cell C contains a manifold with equation 



Z{0 = (7rt)-"V^exp (e - r/^) - /4t 




= a 
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for a large enough value of a. (This property remains valid when the decoherence coefficients 
g^^ are not constant but smoothly varying). 

This means that the normalized mixed states with density matrix P/TrP satisfy the 
criteria for einselected states. This includes their sifting through predictability sieves [TUj . 
since Tr'E?{t) ^ 1 for a reduced density matrix S such that S(0) = P/TrP. The sketch of 
the proof consists in separating diagonalization and mixing according to Eq. (6.9) through a 
canonical transformation maximizing the rate of diagonalization. The sifting property follows 
when t is such that diagonalization has already taken place in the cell C whereas mixing has 
not spilled outside C. This is valid for non-degenerate and degenerate decoherence. 

One can then identify einselected states in general with the classically meaningful states, 
which are defined either as classical properties through the projections P [Hnilll] or as quan- 
tum states by the density operators P/TrP. The predictability sieve criterion is universally 
valid. Its stability under a change of definition for the collective subsystem and the environ- 
ment (i. e. under a shift of Heisenberg's frontier) cannot be proved however along the same 
lines as long as no objective definition of the collective observables has been found. 

9 Conclusions 

As suggested by its name, decoherence is a loss of correlation between local phases of 
a system involving a large number of constituents. It may take in principle many different 
aspects because "in principle" the set of states of a quantum system is enormous, even much 
more than its Hilbert space. Empirical physics is however interested in the systems really 
occurring in nature or built in the laboratory, which can be measured or observed. 

A wide gap between theory and practice is our unability to characterize mathematically 
these "real" systems JIJ'. There is a wide agreement that they always involve some "col- 
lective" degrees of freedom but the problem of their definition from first principles is not 
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yet solved. The study of decoherence will probably remain semi-empirical as long as the 
program suggested in Section 8, or an equivalent one, is not completed. 

The practical results of the present study were concerned with the three main aspects 
of decoherence: suppression of macroscopic interferences, einselection and later classical 
behavior. 

The suppression of macroscopic interferences is a general feature. The interference terms 
disappear for two collective wave functions with a large enough difference in the average 
values of position or momentum (or both). 

Einselection is the election of definite states representing exclusive events with well- 
defined probabilities. It is essential in measurement theory and its properties were the main 
purpose of this paper. Two different cases had to be distinguished, which were respectively 
called degenerate and non-degenerate. 

There is something puzzling in this dichotomy if one does not distinguish also between 
what is most general (or frequent) either from the standpoint of a mathematical theory or 
of empirical physics. A very large class of physical systems leads to the degenerate form of 
einselection, which is practically a diagonalization of the reduced density matrix in the basis 
originating from the collective position coordinates. These systems are truthfully described 
by hydrodynamical variables after coarse graining I13j. Although this condition is still 
restrictive, it turns out in practice that the mechanical parts of a physical system, which are 
described by these variables are entangled most often with other degrees of freedom so that 
degeneracy (with diagonalization) is extended to them. A simpler way of saying this is that 
most observations and measurements involve or could involve a reading of the position of 
some mechanical "pointer" , imposing diagonalization as the outcome of decoherence. 

A sufficient condition for degeneracy is given by the well-known Eq. ()1.3|1 . which covers 
the hydrodynamical case. It is very restrictive however from a mathematical standpoint 
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and, in the absence of a criterion defining a realistic system, one had also to investigate 
the general case of non-degenerate decoherence. The results did not quite agree with the 
conjecture of predictability sieves [10^. One found a tendency of decoherence to combine 
an approximate diagonalization with a partial lumping of probabilities rather than a clear 
mutual exclusion of events, which is typical of unprecise measurements. 

The situation was clearer when one looked at the classical behavior after decoherence. 
There is a simple way to reconcile the present results with the Zurek's predictability criteria 
pnj . It consists in identifying the einselected states with the mixed states representing 
classical properties, which I proposed earlier The fact that these states are best defined 
by the mathematics of microlocal analysis [25j, as well as Fefferman's promising approach 
to the definition of collective observables jSZj , indicate in my opinion that this framework is 
the right one. 
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A - Appendix A. Decoherence as an irreversible 
process 

Auxiliary densities 

The simplest way for computing the quantities Sj for decoherence consists in using the 
observables {\k >< k'\ + \k' >< k\) (g) Jg and {l/i){\k >< k'\ - \k' >< k\) (g) Je (for k ^ k', 
k > k') together with the diagonal terms \k >< k\ ® I^.. They will be denoted altogether by 
A' and they satisfy the relations 

Tr.A^A^' = 5^^' (j,jVl,e). (A.l) 

Rather than, the exponential form ()2.4|1 for the test density operator, it is convenient to 
write it as 

Po = ( E ^'1 ®Pe- (A.2) 

From Eqs. ()A.1|) . one sees that the coefficients in Eq. ()A.2|) are the average values of the 
observables AK Therefore 

dp^/da^ =A^®p, (jVl,e). (A.3) 

After writing po = exp(— a — 13 He), one obtains 

dp /BE = -A^ , da/dE = E/A^ , dpe/da = -p^ , dpjd(3 = -H^Pe (A.4) 

(where A is defined by Eq. fl3.8j) ). In view of the definition ()2.5|) for the auxiliary densities 
(or equivalently the definition ()2.8|) of the projection P), their expression ()3. 5113 .91) in Sec- 
tion 3 follows immediately from Eqs. ()A.3IIA.4]) . 
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The evolution equations 



In the first evolution equation ()2.1H|1 . one can compute PLpi by applying the expression 
()3.9|) giving P to Lpi = —{i/h)[H, pi]. Two traces tr[H,pi] and Tr{He[H, pi]) enter in the 
result. Using cyclic invariance of traces together with Eq. ()3.1U|) specifying Hi and 

Pl= pr^Pe , (A.5) 

one finds that 

tr [H, pi] = [H^, p,] , 

Tr {H, [H, pi]) = Tr {[H,, H] , pi) = Tr {[H,, H^] , pi) . 

Therefore 

PLpi = {-i/h) pr] ®pe + lc® Pe ^"~^ ^ Tr {H, [Hi, pi])| . (A.6) 

In order to compute Ppi, one remarks that although the quantities pe, E and A in Eq. 
()3.9|) are time-dependent, the quantity Tr{HePi) — ETrpi vanishes so that one has 

^ 

Ppl = Pr®Pe- Eh ® Pe ■ (A.7) 

An identical result is obtained for Pp, so that 

Pp2 = . (A.8) 
One thus get the first evolution equation 

pi = PLpi + Ppi + PLp2 , (A.9) 
with the expressions ()A.6|) and ()A.7|) for the first two terms whereas PLp2 is given by Eq. 
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It is convenient to split this equation into one for and another for pe (or equivalently 
for E). This is obtained by taking respectively the trace of Eq. ()A.9|) with respect to the 
enviro nment and the collective Hilbert space. The environment trace of the second term in 
the right-hand side of Eq. ()A.6|) vanishes, as well as trPpi (because trpe = 0). According 
to Eq. ()3.9|) . the environment trace of PLp2 reduces to 

-^tr[H,p2] = -^tr[H,,p2] , 

where the second equality results from the vanishing of tr[Hf., P2] (as the trace of a commu- 
tator) and of tr[Hc, P2] (because of trp2 = 0). One obtains thus the basic equation 

Pr-^{[H,,pr]+tr[H,,p2]) . (A. 10) 

The trace of Eq. ()A.9|) on the collective Hilbert space reduces to a (potentially infinite) term 
{TrcIc)Pe{He — E), multiplied by a number, which must then vanish so that 

E = ~Tr{H,[H^,p, + p2]) = -^Tr{H,[H,,p2]) , (A.ll) 

the last equality resulting from Tr{He[Hi, pi] = Tr{Hi[He, pi], whereas 

[He,pi\=0, (A. 12) 

since pe is a function of Hg. 

The second evolution equation 

One can now write down the evolution equation ()2.14|) for p2, which is formally. 

P2 = QLp2 + QLpi- Ppi (A. 13) 

(after taking Eq. ()A.8|1 into account). This will be done according to Section 4 by considering 
Hi as a perturbation. One needs only to compute p2 at first order in Hi and some terms in 
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Eq. ()A.13|) can be therefore immediately neglected. For instance -Ppi, as given byEq. (1X71) 
is negligible because E is of second order (according to Eq. ()A.11|) and furthermore, 

Pe = -$ {H, -E)p, = EA~^ {H, - E) (A. 14) 

is also of second order (the second equality resulting from Eq. ()A.4|) . One can also neglect 
the commutator [Hi, P2] in Lp2 as being of second order. 

Let us now consider the quantity QLpi = Lpi — PLpi. One has 

Lpi = -{i/h) [H,pi] = -{i/h) [H, + Hi,pi] , (A.15) 

where the second equality results from Eq. ()A.12j) . Then 

PLpi = tr {Lpi) ® pe + (/^ (g)He-E) pi {tr (/^ ® HeLpi) - ETr (Lpi)} . (A.16) 

The last term Tr{Lpi) vanishes (as a trace of a commutator). The preceding term also 
vanishes since 

Tr (/, (g) H.Lpi) = ~{i/h)Tr (/, ® [H, pi]) = {i/h)Tr {H (/, ® H^, pi]) 

and the last commutator vanishes. Therefore 

PLpi = {-i/h)tr [H, pi] 0pe = i-i/h) [JHc, Pr] ®pe = i-i/fi) [He ® /e, Pi] , 

where the second equality results from Eqs. ()3.10|) and ()A.12|) . The first term in Lpi as 
given by Eq. ()A.15j) is therefore cancelled and one is left with the simple equation 

P2 = -{^/h) [Ho, P2] - {^/h) [Hi, pi] (A.17) 

where Ho is the uncoupled hamiltonian: 

Ho = He0le + Ic^He . (A.18) 
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B - Appendix B. Decoherence and dissipation 

Derivation of Eq. I[5.^^) 



According to Eq. (j5.2|) . one must evaluate at lowest order in h the symbol D of the 
collective operator 

D = -{l/h)Hr[Hi,[H^,U{t-t')pr{t')0 Pe{t')U~\t-t')] . (B.l) 

Since U{t) = Uc{t) ®Uf,{t) (with Uc{t) = exp{—iHct/h) and Ue{t) = exp{—iHf,t/h)), one can 
slightly simplify the density operator by writing 

U{t - t')pr{t') ® Pe{t')U-\t - t') = U,{t - t')pr{t')U;\t - t') ® Pe(t') , (B.2) 

in view of the equality (resulting from Eq. ()A.14jl and valid up to order Hi) 

U,{t-t')p,{t')U~\t-t')^p,{t') . 

Letting A and B be two arbitrary operators, A{x,p) and B{x,p) their operator-valued 
symbols, Eq. ()5.5|) gives the symbol of the commutator [A, B], with the notation of Section 5: 

{-th/2) (ApB, - A^Hp - BpA, + B,Ap) + Oih') , (B.3) 

Eq. ()B.3|) can be used twice for obtaining the symbol of the double commutator in Eq. 
dHU). The symbol of the operator Uc{t - t')pr{t')U-^{t - t'), which will be denoted by W'^ 
is an ordinary function and it commutes with the operator- valued symbols Hi , H^ and their 
derivatives. After a straightforward calculation, one gets: 

D = J' dt' (^[C'^'^Wj + C''''W^) +^[CP''Wj + CPPW^)^ . (B.4) 
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The coefficients are given by 

C^^ = Ur { (h,hI + hIh,) p,} , (B.5) 

= -\tr { (pr^l + hIh,) pe} , (B.6) 

C'^' = —tr { {hJiI + hIh.) Pe} , (B.7) 

= Ut { {hJiI + hIh:) p,} , (B.8) 

Explicit expressions 

The decoherence coefficients after neglecting retardation are obtained by integrating the 
coefficients ()B.5IIB.8|) on the time t' . It is convenient to introduce the matrix elements of 
Hi : 

<n\H^^\n' >= Hunn' exp {-iunn'r) 
where the index i denotes either x or p, t = t — t' and Hunn' =< n\Hii\n' >. One has then 

^tr I (= HiH^ + H^jH^ Pe} = XI Hiinn'Hijn'n exp {iuJnn'r) Pnn' COsh {(3hUnn'/'2) ■ 

nn' 

The decoherence coefficients 

Jo 

are then given by 

roo 

= rfr X Hipnn'Hipn'n exp (iuJnn'T) Pnn' COsh (PkUnn' f^) , (B.9) 

^ nn' 

POO 

gPX ^ gxp ^ _ Hipnn'Hi^n'n eXp {iuJnn'T) Pnn' COsh {pflUnn'/'^) , (B.IO) 

JO I 

nn' 

roo 

gPP = / dr V-ffiw-f^ixn'nexp (zu;„„/r)p„„/ cosh(/5/iu;„„'/2) • (B-H) 
Jo I 

nn' 

The quadratic form in two real variables (a, a') with these coefficients is given by 

.g^'^C? + 2g'^^0L0i + g^'^ol'^ = drY^ IHi^nn'a - Hip^'n^' cos {uJnn'T) Pnn' COsh {PhuJnn'/2) 

^ / 

(B.12) 
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nn' 



and it is clearly non-negative. 



Dissipation coefficients 



Let now A^. {k = 1 or 2) denote either X or P. The time derivatives of their average 
values are given by 

^^A^ = ^Tr i[H, A,] p) = ^Tr^m, A,] p,) + '-Tr i[H^, A,] p,) . (B.13) 

The first term is due to collective dynamics and the second one represents dissipation effects. 
Using Eq. ()4.5|) for this dissipative term reads 

- ^ drTr {[ifi, Au] [H^ , U{T)pr ® PeU'\r)\ } . (B.14) 

Eq. ()5.7p can be used to replace the collective trace by an integration over phase space of 
the operator symbol. One can use 



When computing the symbol of [H'[,U{T)pr ® PeU^^i^r)], one will retain only the term 
originating from the Poisson bracket between Hf and either U{t) or U^^{t), because it can 
be seen that all the other contributions do not contribute to the phase space integral or are 
of higher order in h. Finally, denoting by Hc{x,p) the symbol of (which is the collective 
Hamilton function), one gets 

[HT,UiT)pr^pM-'{r)] = y { {Hlpe - PeHl) V'ix) - {hIp, - pjll) tW. 

In the semi-classical case, when the Wigner function is slowly varying, Eq. ()B.13|1 generates 
the classical equations of motion, which read (after writing {p,x) in place of < P >, < X > 
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and neglecting retardation): 

^ = -H,^ - Y^Hcp - Y'^Hc. , (B.15) 

(j'T 

— = -H,p - j-^H^p - . (B.16) 

The coefficients are explicitly given (after performing a partial integration over the time 
r) by 

roo 

= dr H i:,nn' H i^n' n i'i'^nn' r) Pnn' ismh{pnUJnn' /huJnn') , (B.17) 

Jo 

= - / drHi^nn'Hipn'n exp iiuJnn'r) Pnn' (sinh {(3huJnn'/2) /flUJnn') = 7''^,(B.18) 

Jo 

roo 

7^^= / dTHipnn'Hipn'new{'i^nn'r)Pnn' isi'ah{PhuJnn'/'2) /hUnn') ■ (B.19) 




Comparing these results with Eq. ()B.9IIB.lT|) . a strong formal similarity appears between 
the dissipation and the decoherence coefficients. They are even directly proportional at high 
enough temperature (when T ^ hQ so that cosh(/5^to'„„'/2) ^ 1 and sinh(/9/ico'„„'/2) pa 
Phunn'/'^), namely: 

^ g'^/T . (B.20) 

One may also notice that the dissipation of collective energy dEc/dt is always negative (or 
zero), since the quadratic form with coefficients Y'' is non-negative for the same reason as 
in Eq. ()B.12|) . Finally, the equality Y^ = 7^^ is a special case of the Onsager symmetry 
relations (24). 
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C - Appendix C. Previous models 



Much of our knowledge on decoherence was first obtained from a study of simple models. 
Since the present theory claims a much wider range, it should at least recover these older 
results. This will be the topic of the present Appendix. 

Oscillator bath 

A model of environment consisting of a bath of harmonic oscillators has been much 
investigated. For an atom interacting with radiation in a cavity, the model is exact and its 
results have been successfully compared with experiment ^Ij. The atom is represented in 
that case by a two-state system. When considering a collective system with position X, the 
coupling Hamiltonian is typically written as 



where the sum is over all the oscillators, Aj being a coupling constant, and a] their anni- 
hilation and creation operators. 

It is easy to use this coupling in the formulas of Appendix B and to recover the pre- 
viously known results. The calculations are essentially trivial and need not be given here 
explicitly. Perhaps more interesting is the question of the range of this model. In addition 
to their thorough study of it, Caldeira and Leggett suggested that it should be considered 
as very general They start from the fact that the number of energy eigenstates \n > 
of the environment is extremely large. They introduce formally an oscillator for each such 
state (i.e. the index i will stand for n in Eq. ()C.1|) ). They notice that a state |r;, > is either 
occupied or not, these two possibilities being represented by the ground state and the first 
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excited state of the associated oscillator. Then they argue that higher excited states of the 
oscillators will play no role because of the smallness of their probability of excitation. They 
conclude that an oscillator bath can represent almost any environment. 

The argument is however erroneous, because the coupling resulting from their proposal 
would not have the very simple form of Eq. ()C1|) in general, because transitions i ^ j with 
i ^ j are certainly essential and they do not appear in the Hamiltonian ()C.1|1 . Moreover, 
the transitions i — > i with a change of occupation number are not correctly represented in a 
quantity such as 

tr (HiHiPe) = < n\Hi\n' >< n'\Hi\n > p^n ■ 

nn' 

If the occupied state \n > is considered as the first excited state \i, 1 > of an oscillator, the 
contribution of n' >= \i, 2 > cannot be omitted from the sum (??) if the Hamiltonian ()C.1|) 
is used. The interest of a theory ignoring the constraints of the oscillator model cannot 
therefore be disputed. 

The scattering model 

Joos and Zeh have proposed a beautiful model of decoherence, when an object with posi- 
tion X interacts with a bath of particles P|. Decoherence appeared then as an accumulation 
of scattering effects. The resulting master equation looked very similar to one occurring in 
the oscillator model and this was a very strong hint for some universality in the mathemat- 
ical expression of decoherence. This universal character has been explained in the present 
paper, but something would still be missing if the scattering model were not also derived. 
This derivation is non-trivial so that it will be given explicitly. 
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Let us consider for definiteness a spherical solid object with center X. As explained in 
Section 3, the pressure exerted on it by the outside gas is included in He so that Hi represents 
the effect of random collisions of the outside molecules (or photons) on the sphere. Rather 
than doing a complete calculation, it will be sufficient to look at one term in Eq. ()5.4j) . 
namely (with h = 1) 

Di = - f' dt' tr {HiUcit - t')pr{t')U-\t - t') ® peH^} . (C.2) 

One may consid er first the case of an environment consisting of particles having the 
same momentum, in a pure state \k > : pe = \k >< k\. The corresponding wave function is 
a plane wave exp(ifc.r) and, with this normalization, there is one particle per unit volume. 
Their flux $ is the velocity u of the particle and a sum over one-particle states \k' > amounts 
to an integration over dk/ (27r)^. 

A few remarks are useful: 

1. One can neglect the time evolution factors Uc{t — t'), which are slowly varying. 

2. Introducing momentum eigenstates \p > of the object and the outgoing scattering 
states \p, k, out > associated with the plane waves \p, k >= \p > ®\k >, one can write 

< p2, k2, out\Hi\pi, ki >=< p2, k2\T\pi, ki> 6 {pi + ki - p2 - k2) , (C.3) 

where T is the (off energy-shell) collision matrix for the scattering of a particle on the object 
[45]. 

3. One can use the invar iance of the scattering T-matrix under a change of reference 
frame. For non-relativistic values of p and taking into account the large mass of the object 
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(as compared with the particle mass), a Gahlean transformation with small velocity gives 

< P2, hlTlpi, ki >=< p2 - pi, k2\T\0, ki >< pi, ki\T\p2, ^2 >=< 0, ki\T\p2 - pi, ^2 > • 

(C.4) 

4. In view of the delta function in Eq. ()C.3|) . it is enough to know ki and ^2 for fixing 
P2 — pi so that the right-hand side of Eq. ()(I4j) can be written more simply T{ki,k2). 
Conversely, using Fourier transforms, one can introduce the states \x,k > corresponding to 
a localized object, whereas out states \x, kout > involve the same value of x and a scattered 
particle. One thus gets 

< x' , k20ut\Hi\x, ki >= T {ki, k2) S{x — x') exp {i {ki — k2) x} . (C.5) 

This result has two important consequences: It shows that Hi acts like a function of X, 
although this property shows up only when scattering states are used. The imaginary expo- 
nential in the right-hand side is moreover typical of the method that was used by Joos and 
Zeh. 

5. In view of Eqs. ()5.3|) . and (jC.3jl . one has 

< x, ki\Hi\x\ k2, out >= exp {iuj{t — t')} 5{x — x') exp {—i {ki — ^2) x} -T* {ki, ^2) , (C.6) 
where a; = E{k2) — E{ki). 

6. When throwing out retardation, the integral on t' of the imaginary exponential in 
(IC.6|) gives 'k5{uj) — iP{l/uj), where P stands for a Cauchy principal part. It may be shown 
however that another term in D originating from piHi cancels the principal part and one 
must keep only therefore the delta part. 
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Finally, the matrix element < xlZ^ilx' > can be easily computed if one uses the orthonor- 
mal set of outgoing states {\k' out >} when computing the trace tr. One gets 

<X\Di\x'>= J [dk/ (27:^) TT6{u)\T{k,k')\^ exp {t{k - k'){x - x')} pr{x,x') . (C.7) 

But the quantity dk{27r)^'^6{uj)\T{k, k')\'^ has a very simple interpretation: It coincides with 
the product da^ of the differential cross-section da for the scattering k ^ k' times the flux 
$ of the environment particles having the given momentum k. One can then replace the 
trivial density matrix \k >< k\ by a thermal density and introduce the various different 
particles in the gas, thus obtaining: 

< x\Di\x' >= J ■Kdad^exp{i{k - k'){x - x')} pr{x,x') . (C.8) 

Similar results are obtained for the three other terms in D but it will not be necessary to 
push the calculation further since, from there on, it becomes identical with the one by Joos 
and Zeh. Their method was of course simpler than the present one, as one expects from an 
intuitive approach compared to a technical one. The present calculation shows however how 
universal and versatile the fundamental master equation ()4.(jj) is. 
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